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Abstract 

In this paper we construct modified cosmic Chaplygin gas which has viscosity. We use 
exponential function method to solve non-linear equation and obtain time-dependent 
dark energy density. Then discuss Hubble expansion parameter and scale factor and 
fix them by using observational data. We also investigate stability of this theory. 
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1 Introduction 



It is believed that the most part of Universe filled with dark matter and dark energy. There- 
fore, dark energy and related topics are important subjects to study in theoretical physics 
and cosmology. An important problem is determining nature of dark Universe. It is found 
that the dark matter may be consists of neutrinos [1] axions [2] or WIMPs (weak interactive 
massive particles) [3]. In that case there are several ways to specify the nature of the dark 
Universe. For example studying time-dependent density help to give information about dark 
matter [4] and dark energy [5]. 

After that we need a model to describe dark Universe. In that case there are some phe- 
nomenological and theoretical models which based on discovery of the accelerating expansion 
of the Universe [6, 7]. Some of the famous phenomenological models for dark energy briefiy 
explained below. 

The cosmological constant and its generalizations are the simplest way to modeling the dark 
energy [8]. Another candidates for the dark energy are scalar- field dark energy models. A 
quintessence field [9] is a scalar field with standard kinetic term, which minimally coupled to 
gravity. In that case the action may has a wrong sign kinetic term (minus instead of plus), 
which the scalar field is called phantom or ghost [10]. While there is a quantum instability 
in the phantom models but this model to be consistent with CMB observations [11]. Com- 
bination of the quintessence and the phantom is known as the quintom, which is another 
model for dark energy [12]. 

Extension of kinetic term in Lagrangian yields to a more general frame work on field theoretic 
dark energy, which is called k-essense [13, 14]. A singular limit of k-essense is another model, 
named Cuscuton [15]. This model has an infinite propagating speed for linear perturbations, 
however causality is still valid. 

The most general form for a scalar field with second order equation of motion is the Galileon 
field which could behave as dark energy [16]. Another extension of these models is called 
the ghost condensation, which also solved the quantum instability of phantom dark energy 
[17]. There are also various studies in holographic dark energy models (see Refs. [18-22]). 
However, presence of a scalar field is not only requirement of the transition from a Universe 
filled with matter to an exponentially expanding Universe. The matter components in cos- 
mology are written in terms of fiuids, so most of dark energy models have fiuid description. 
Therefore, Chaplygin gas (CG) used as an exotic type of fiuid, which is a model for dark 
energy [23, 24]. This model based on Chaplygin equation of state [25] to describe the lifting 
force on a wing of an air plane in aerodynamics. The CG was not consistent with observa- 
tional data of SNIa, BAO, CMB, and so on [26-29]. Therefore, an extension of CG model 
proposed [30, 31], which is called generalized Chaplygin gas (GCG), and indeed proposed 
unification of dark matter and dark energy. However, observational data ruled out such a 
proposal. Then, GCG extend to the modified Chaplygin gas (MCG) [32]. There is still more 
extension such as generalized cosmic Chaplygin gas (GCCG) [33]. In this paper we deal with 
next extension which is modified cosmic Chaplygin gas [34]. 

On the other hand bulk viscosity plays an important role in the evolution of the Universe. 
The idea that Chaplygin gas may has viscosity first proposed by the Ref. [35] and then 
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developed by [34, 36, 37, 38, 39]. 

In the Ref. [35], for the first time, GCG with bulk viscosity studied. In the Ref. [34] we in- 
deed extend Ref. [35] to the case of viscous MCG for the special case of a = 1/2 in equation 
of state formula and /c = in FRW metric. Then, in the Ref. [37] we extend our previous 
work to the case of MCCG and calculate time-dependent dark energy density and discussed 
stability of system. In the Ref. [38] viscous CG in non-fiat FRW Universe has been studied 
and time-dependent density of dark energy obtained under assumption of cosmic expansion 
without acceleration or deceleration which means power law scale factor. Then in the Ref. 
[39] we extend our previous work to the case of arbitrary a instead of a = 1/2 and studied 
viscous GCG. In all cases the time- dependent density calculated approximately, and under 
some assumption for simplicity. 

Now, in this paper we would like to study viscous MCCG with arbitrary a and use powerful 
tools to obtain more exact solutions. Indeed, in order to solve non-linear equation we use ex- 
ponential function method which is more appropriate than previous approximate methods. 
Furthermore we investigate evolution of scale factor and try to compare our results with 
observational data, and also discuss stability of system by using speed of sound in fluids. 
This paper is organizing as the following. 



2 Review of Chaplygin gas model 

One of the recent cosmological models which is based on the use of exotic type of perfect 
fluid suggests that our Universe filled with the Chaplygin gas with the following equation of 
state [38, 40], 

P = --, (1) 

where i? is a positive constant. The equation (1) introduced by Chaplygin to describe the 
lifting force on an airplane wing [25]. Chaplygin gas is also interesting subject of holography 
[41], string theory [42], and supersymmetry [43]. It is also possible to study FRW cosmology 
of a Universe filled with CG [44]. CG equation of state (1) has been generalized to the form 
[24, 45, 46, 47], 

V= , 2 

with < a < 1, which is called GCG, and is also interesting from holography point of view 
[48]. As we can see from the equation (2), the GCG is corresponding to almost dust (p = 0) 
at high density which is not agree completely with our Universe. Therefore, MCG with the 
following equation of state introduced [32, 49, 50], 

where 7 is a positive constant. This model is more appropriate choice to have constant 
negative pressure at low energy density and high pressure at high energy density. The 
special case of 7 = 1/3 is the best fitted value to describe evolution of the Universe from 
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radiation regime to the A-cold dark matter regime. Also, the MCG with 7 = 0.085 and 
a = 1.724 is coincident with some observational data [51]. The next extension performed by 
the Ref. [33] where the GCCG introduced by the following equation of state, 



(4) 



where C is a constant. In the case of w = one can write C = B — 1. The speciality of 
this model is stability so the theory is free from unphysical behaviors even when the vacuum 
fluid satisfies the phantom energy condition. In this paper we construct the MCCG with the 
following equation of state. 



P 



IP 



B 



ip 



l+a 



B 



(5) 



Our main goal is using above equation of state to study FRW bulk viscous cosmology, 
therefore we need to review FRW bulk viscous cosmology and then add viscosity to MCCG. 



3 FRW bulk viscous cosmology 

As we know the Friedmann-Robertson- Walker (FRW) Universe is described by the following 
metric, 

ds' = + a(t)2(-^^ + r'dn'), (6) 
1 — kr^ 

where dVt^ = d9^ + sin^ Odcj)"^, and a{t) represents the scale factor. The 9 and parameters 
are the usual azimuthal and polar angles of spherical coordinates, with < < tt and 
< < 2tt. The coordinates {t, r, 6, 0) are called co-moving coordinates. The constant 
k defined curvature of space so, k = 0,1, and —1 represents fiat, closed and open spaces 
respectively. Our interest in this paper is the first case namely k = 0. In that case the 
Einstein equation is given by, 

Rfiu '^dfj.i'R T^i,, (7) 

where we assumed c = 1, SvrG = 1 and A = 0. It is assumed that our Universe is filled 
with the MCCG which plays role of dark energy with equation of state (5). Using the line 
element (6) and the Einstein equation (7), the energy- momentum tensor corresponding to 
the bulk viscous fluid is given by the following relation, 

T^u = (P + p)u^Uy - pg^y, (8) 

where p is the dark energy density and is the velocity vector with normalization condition 
u^Ui, = —1. Also, 

p = p-3CH, (9) 
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is the total pressure which involves the proper pressure p, given by equation of state (5), 
bulk viscosity coefficient ( and Hubble expansion parameter H = a/a. It is also assumed 
that the dark energy is conserved with the following conservation equation, 

p + 3H{p + p) = 0. (10) 

In the next section we write field equations of viscous MCCG, and then try to solve them 
for dark energy density. 



4 Viscous modified cosmic Chaplygin gas 

By using the relations of previous section one can obtain the following field equations 

P 
3' 

and 



= fill 



H + H' = -^-^, (12) 
6 2' ^ ' 

where dot denotes derivative with respect to cosmic time t. The energy-momentum conser- 
vation law (10) reduced to the following equation, 

p + 3{p + p)H-9CH^ = 0, (13) 

where d^T^'^ = is used. The equation (11) can be easily extended to the case of non-flat 
Universe and non- vanishing cosmo logical constant, 

p k A 
3 ~ 3" 

As we can see the second term of right hand side contains scale factor which complicate our 
calculations, therefore in this work we consider the simplest case with k = and A = 0. 
Now, we would like to solve equation (13) to obtain time-dependent energy density. Using the 
equation of state (5) in the energy-momentum conservation formula (13) gives the following 
differential equation. 



a 

p + 3- 
a 



^'^^^P-JU^ 7 '^a 



0. (15) 



Then, by using the field equation (11) and change of variable p = X"^, one can rewrite the 
equation (15) approximately as the following non-linear differential equation, 

X + biX^ + + 63X = 0, (16) 

where dot denotes derivative with respect to time, and we defined, 

/3 

bi = — (1 + 7 + c^), 

h = -^C- (17) 
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5 Time-dependent energy density from exponential func- 
tion method 



According to the exponential function method [52], the nonhnear differential equation (16) 
has the following solution, 

x = 'j^:l±^i±^^ (18) 

where coefficients q can be obtained by substituting the solution (18) in the equation (16). 
This gives the following coefficients, 

_ (2&i6 - 2hhh -hl + l)c + (4&i6 + 2b3 - 2)e""* 

26i(2M + &3 + l) ' ^ ' 

(Abjb% + ^hhhl + 4&3& + ^hlh^ + hl + hl-h- l)c 
~ 261(2616 + 63 + 1)2 

(46^636^ + 4616^6 + hl + 4blb^ - 46i6 - 6^ - 63 + l)e-"* 
+ 261(2616 + 63 + 1)2 ^ ' ^ ^ 

_ (1 - 461663 - 62 - 46262)c + (46i + 263 - 46^62 - 463616 - 6^ - l)e""* 

(2616 + 63 + 1)2 ' ' ^^^^ 

with, 

6=(i^)- 



(23) 



also c and k are free parameters. If c = then any dependence on time vanishes and dark 
energy density reduced to a constant and will be conserved during time. So, in order to 
obtain time- dependent dark energy density we consider the case of c = 1. On the other 
hand expanding Universe suggest that k = 1 to have infinitesimal density at t — )■ 00 limits. 



Hence, one can write. 



c.e- + + cse 
^ [ e-* + l + C4e* J ' ^ ' 

In the Fig. 1 we draw typical behavior of energy density in terms of time for various models. 
Yellow line is corresponding to VGCG which is lower than other cases. Extension of this 
model to VMCG denoted by cyan line. We see that the effect of 7 is increasing energy 
density. VGCCG presented by Red line. We take values of constants 7 and B from the Ref. 
[50]. Finally our interesting case of VMCCG denoted by green line which has bigger energy 
density of other cases. 

We find that the effect of viscosity is increasing of dark energy density. On the other hand 
decreasing of parameter a increases the value of dark energy density. It is also found that 
increasing of cosmic parameter u increases the value of the dark energy density. 
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In all expected, the dark energy density is decreasing function of time which is 

agree with expansion of Universe. 

In summary, we can conclude that the evolution of the dark energy density with time is 
faster than other models. 

In the Fig. 2 we give plot of the dark energy density in VMCCG model for a = 0.1 which 
suggests that the value of the dark energy density is very high. 



p 




Figure 1: Time-dependent dark energy density. Green: VMCCG with 7 = 0.3, and B = 3.4 
[50], u = 0.5, a = 0.5 and C = 1- Red: VGCCG with 7 = 0, and 5 = 3.4 [50], u = 0.5, 
a = 0.5 and C = 1- Cyan: VMCG with 7 = 0.3, and B = 3.4 [50], w = 0, a = 0.5 and C = 1- 
Yellow: VGCG with 7 = 0, and = 3.4 [50], = 0, a = 0.5 and C = 1- 
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Figure 2: Dark energy density of VMCCG with 7 = 0.3, and B = 3.4 [50], u = 0.5, a = 0.1 
and C = 1- 



6 Hubble parameter and Scale factor 

By using the equation (24) and (11) we can find behavior of the Hubble expansion parameter. 
Numerically we draw Hubble parameter in terms of time in the Fig. 3 (left), and in terms 
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of viscous parameter ( in the Fig. 3 (right). We can see that there is a singular point for a 
suitable value of viscous coefficient. 

Time-dependent energy density (24) suggests that the scale factor takes the following general 
form, 

a{t) = (me* + n)'^e-"*, (25) 

where we used the relation H = a/ a. Equation (25) contains new constant parameters as 
m, n, (5, and a which will be fixed by using observational data. We draw typical behavior of 
scale factor (25) in terms of time in the Fig. 4. We can see that Hubble expansion parameter 
and scale factor are decreasing function of time. 

Before we try to fix new parameters and obtain scale factor exactly, it is useful to check 
stability of this theory which is subject of the next section. 




10 0.5 0.6 0.7 0.8 0.9 1.0 

t c 



Figure 3: Hubble expansion parameter in VMCCG model with 7 = 0.3, B = 3.4, 00 = 0.5 
and a = 0.5. Left: In terms of time with C = 1- Right: In terms of ( with t = (Red line), 
t = 1 (Green line), and t = 5 (Blue line). 




Figure 4: Scale factor in VMCCG model. 
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7 Stability 

There are several ways to investigate stability of a theory. In this paper we use reality of 
sound speed, 

C's=f> 0. (26) 
dp ~ ^ ^ 

Numerically, we draw in terms of time in the Fig. 5. It shows that the VMCCG model 
of cosmology is completely stable and there is no instable region. In the Fig. 6, we draw 
in terms of u and (. The left plot of the Fig. 6, shows that increasing u reduced the sound 
speed. 




Figure 5: Square of sound speed in terms of time. 




Figure 6: Square of sound speed with 7 = 0.3, B = 3.4, and a = 0.5. Left: In terms of 00 
with C = 1- Right: In terms of ( with u = 0.5. t = 0, t = 1, and t = 5 represented by Red, 
Green and Blue lines respectively. 
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8 Observational data 



First of all we consider current value of Hubble expansion parameter which is Hq 68 — 72 
[53]. Therefore, according to the right plot of the Fig. 3 (red line), visous coefficient ( should 
be about 0.6 or 0.75, so our results agree with observational data. 

In order to fix the scale factor (25) we use the following observational parameters [54]. 
Declaration parameter is given by the following relation, 

^ (27) 



d2 dt^ 

and jerk parameter is given by, 



then, snap parameter obtained by using the following relation. 



also, lerk parameter is given by. 



s = --— (29) 

d^a , , 

^ = TT^. 30 

a5 dt^ ^ ' 

Current value of these parameters Qq, Jq, Sq, and /q appear in the following Taylor expansion 
around ao, 

^ = 1 + iJot - IqoH'ot' + ^joH'.t' - j^soH^t' + ^kHlt' + ■■■. (31) 



In that case on can obtain 



mnS 

go = -1 - ^ ^- (32) 

(m(o — a) — any 



Following we use results of the Refs. [54-58] to fix our parameters. 



8.1 SNela 

According to the SNela observational data [55, 56] the current value of the declaration 
parameter may be go = ~1- This situation obtained by three separated ways of m = 0, 
n = or 5 = 0. In all cases one can obtain q = s = 1 and j = I = —1. It means that 
observational parameter have no any time dependence. It should be noted that = —1 and 
Jo = 1 are completely agree with the Refs. [55, 56]. 
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8.2 ACDM 

According to the ACDM observational data [57, 58] the current values of the declaration 
and jerk parameters may be go = —0.6 and jo = —1 respectively. Scale factor (25) in the 
jerk parameter (28) gives the following expression, 

mn6{3na — 3m{6 — a) ~\- m — n) 

Jo = 1 7 7T ^ 7^ • (33) 

(m(o — a) — any 

Assuming a = 6 gives us, 

go = -1 - 34 
n 

which together with go = —0.6 yields to, 

m , , 

— = -0.4cT. 35 
n 

Combining the relation (35) with the current jerk parameter (33) together with observational 
data jo = ~1 tell us that o ^ —0.4. Therefore on can obtain, 

0.4t 

a(t) = ^. (36) 

^ ^ [n(0.16e* + 1)]°-' ^ ' 

With n = 0.1 the current value of scale factor coincide with the Fig. 4. 



8.3 Best fitted 

According to the best fitted parameters of the Ref. [54] we know go = —0.64 and jo = 
1.02. These observational data together with the equations (32) and (33) give the following 
equation, 

m^ + n^ + (255 - 2)mn = 0, (37) 
which suggests the following solution, 

m _ 255 - 2 
'n ~ 2 

We can use this relation with variation of 6 to obtain more exact solution. In the table 1, 
we can see that 6 = 2 gives the best solution relative to the best fitted data [54], and yields 
to obtain the following time-dependent scale factor, 

a{t) = a(0)(l - 0.02e*)2e-°-2^^*, (39) 

where we assumed a(0) = n^. 



-1 ± 



4 



(25(5-2)2 



(38) 
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Table 1: Best fitted values of the deceleration and jerk parameters [54], and our theoretical 
values. It shows that 5 = 2, and a = 0.299 have most agreement with observational data. 



9 Conclusion 

In this paper we constructed viscous modified cosmic Chaplygin gas. We reviewed evolution 
of Chaplygin gas to modified cosmic Chaplygoin gas which considered as dark energy, and 
written corresponding equation of state. Then, we introduced viscosity to the system and 
obtained modified conservation equation and field equations. These yield to a non-linear 
differential equation which gives the time- dependent dark energy density, we solved this 
equation by using exponential function method, and study behavior of the dark energy 
density. Then we calculated scale factor and discussed Hubble parameter numerically. We 
also investigated stability of the theory and found that viscous modified cosmic Chaplygin gas 
model is completely stable and there is no instable region. Calculating the scale factor yields 
to unknown coefficients which determined by using the observational data of deceleration 
and jerk parameters. In order to find more exact solutions one can consider also snap and 
lerk parameters. 

In this paper we considered the case of fiat space {k = 0) and zero-cosmological constant. 
It is possible to include these parameters. In that case the field equation modified as the 
following, 

H^ = P-^ + ^. (40) 
3 a2 3 ^ ^ 

In order to study effect of space curvature and cosmological constant one should use the 

equation (40) in the relation (13) to obtain dark energy density. The most problem is the 

second term of right hand side. So, the case of k = and A 7^ is straightforward. Just 

enough to replace p with p + A. 

Acknowledgments It is pleasure to thank M. R. Ganji and A. A. Seifi for introducing 
methods of solving non-linear differential equations. 
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